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Abstract 

The  growing  interest  in  giant  magnetostrictive  materials  for  generation  of  strains  in  smart  structure  systems 
motivates  the  development  of  increasingly  accurate  performance  prediction  and  optimization  tools.  We  propose  a 
model  for  the  strains  generated  by  magnetostrictive  materials  in  response  to  applied  magnetic  fields.  The  direct  or 
magnetostrictive  effect  is  modeled  by  considering  active  and  passive  components  of  the  strain.  The  active  or  external 
component,  associated  with  the  alignment  of  magnetic  moments  with  the  external  magnetic  field,  is  modeled  with 
a  ferromagnetic  hysteresis  model  in  combination  with  a  quartic  magnetostriction  law.  The  passive  or  internal 
component,  associated  with  the  elastic  response  of  the  transducer  materials  as  they  vibrate,  is  modeled  through  force 
balancing  which  yields  a  wave  equation  with  magnetostrictive  inputs.  The  effect  of  stress  on  the  magnetization  of 
the  magnetostrictive  core,  or  the  magnetomechanical  effect,  is  implemented  by  considering  a  ‘law  of  approach’  to  the 
anhysteretic  magnetization  caused  by  stress.  This  provides  a  representation  of  the  bi-directional  coupling  between 
the  magnetic  and  elastic  states.  It  is  demonstrated  that  the  model  accurately  characterizes  the  magnetic  hysteresis 
and  the  strains  output  by  a  prototypical  transducer. 
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1  Introduction 

This  paper  presents  recent  advances  on  the  modeling  of  strains  generated  by  transducer  systems  which  utilize 
magnetostrictive  materials  as  the  active  element.  By  virtue  of  the  reciprocal  nature  of  magnetostrictive  materials, 
such  transducer  systems  are  capable  of  providing  both  actuation,  characterized  by  magnetic-to-elastic  energy  trans¬ 
duction,  and  sensing,  characterized  by  elastic-to-magnetic  energy  transduction.  One  crucial  aspect  of  this  two-way 
energy  exchange  is  that  neither  transduction  mechanism  occurs  independently,  but  rather  they  both  occur  simul¬ 
taneously  during  transducer  operation.  In  the  direct  or  magnetostrictive  effect  (James  Joule,  1842),  the  action  of 
a  magnetic  field  and  ensuing  magnetization  induces  strains  in  the  magnetostrictive  core.  These  strains  are  in  turn 
associated  with  a  stress  field  which  affects  the  magnetic  state,  effect  known  as  the  inverse  magnetostrictive  effect  or 
simply  the  magnetomechanical  effect. 

The  close  connection  between  magnetostriction  and  the  magnetic  behavior  under  stress  has  been  long  recognized, 
and  extensive  experimental  evidence  on  the  magnetomechanical  effect  has  been  documented  in  the  literature  [1,  2,  3]. 
In  recent  years  there  has  been  renewed  interest  on  this  phenomenon  [4,  5]  because  of  its  relevance  for  technological 
applications  such  as  non-destructive  evaluation  and  sensing.  Interest  in  the  effects  of  stress  on  the  performance  of 
highly  magnetostrictive  materials,  such  as  Terfenol-D,  is  also  motivated  by  the  need  for  increased  levels  of  adaptability 
in  smart  structure  systems. 

In  this  paper,  we  extend  the  model  derived  in  [6,  7],  which  includes  magnetic  hysteresis  and  nonlinearities  as 
well  as  structural  vibrations,  to  incorporate  the  magnetomechanical  effect  in  a  magnetostrictive  rod  as  it  drives  a 
transducer.  The  model  is  illustrated  in  the  context  of  the  Terfenol-D  transducer  depicted  in  Figure  1,  but  is  not 
limited  to  this  application.  The  present  approach  may  be  implemented  for  other  transducer  materials  such  as  nickel. 
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Because  the  model  quantifies  the  full  magnetomechanical  coupling  exhibited  by  the  material,  it  provides  a  more 
comprehensive  representation  of  the  physics  of  the  system  than  does  the  uncoupled  model.  The  model  is  derived 
under  the  assumption  that  the  transducer  is  operated  in  quasistatic,  temperature  independent  conditions  and  that 
the  material  is  maintained  in  an  approximately  isotropic  state. 

The  uncoupled  formulation,  modified  to  accommodate  additional  terms  arising  from  the  overall  coupled  approach, 
is  summarized  in  Section  2.  The  magnetomechanical  effect  model  employed  is  discussed  in  Section  3,  while  the 
complete,  coupled  model  is  compiled  in  Section  4.  In  Section  5,  we  compare  model  results  with  experimental  data. 
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Figure  1:  Cross  section  of  the  prototypical  Terfenol-D  transducer  employed  in  the  model. 


2  Magnet ostrictive  Effect 

In  previous  investigations,  the  magnetostrictive  effect  has  been  modeled  in  three  stages  [6,  7].  In  the  first  stage, 
application  of  a  magnetic  field  H  produces  a  bulk  magnetization  M .  The  hysteresis  and  saturation  effects  present 
in  the  M  —  H  relationship,  illustrated  in  Figure  2(a),  are  modeled  through  a  mean  field  model  of  ferromagnetic 
hysteresis  originally  developed  by  Jiles  and  Atherton  [8,  9,  10].  The  second  stage  involves  the  characterization 
of  the  magnetostriction  due  to  the  alignment  of  magnetic  moments  with  the  applied  magnetic  field  and  ensuing 
magnetization.  We  call  this  magnetostriction  A,  to  distinguish  it  from  the  total  magnetoelastic  strain  e.  The  most 
salient  characteristics  of  the  A  —  M  relation  are  depicted  in  Figure  2(b).  The  model  employed  at  this  stage  is  a 
truncated  empirical  series  expansion  symmetric  about  M  =  0.  The  third  aspect  that  needs  to  be  characterized  is 
the  total  strain  e,  which  incorporates  both  magnetic  effects  and  the  elastic  properties  of  the  transducer.  Typical 
issues  to  address  in  the  £  —  H  relationship  are  hysteresis,  nonlinearities,  and  saturation  effects,  as  indicated  in 
Figure  2(c).  Modeling  of  the  structural  vibrations  of  the  magnetostrictive  core  is  posed  through  force  balancing,  in 
the  form  of  a  PDE  equation  which  includes  the  intrinsic  magnetostriction,  system  compliance,  internal  damping,  and 
corresponding  boundary  conditions  associated  with  the  mechanical  transducer  design. 

It  is  noted  that  the  rate  of  change  of  magnetization  with  time  can  be  expressed  as 

ciM  _  f  8M  \  dH  ( dM  A  da 
dt  \  dH  )  dt  \  da  )  dt  ’ 

which  motivates  treating  the  field  and  stress  components  independently.  The  remainder  of  this  section  deals  with 
the  field-dependent  magnetization  and  corresponding  derivative 
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(a)  (b)  (c) 

Figure  2:  Model  simulations  depicting  stages  1-3  of  the  magnetostrictive  effect:  (a)  constant-stress,  relative  magnetization  M- 
versus  applied  field  H,  with  Me  the  saturation  magnetization,  (b)  relative  magnetostriction  versus  relative  magnetization, 
with  As  the  saturation  magnetostriction,  and  (c)  total  strain  e  versus  applied  held  H. 


2.1  Field-Dependent  Magnetization 

As  discussed  in  [11,  12],  the  magnetization  process  taking  place  in  ferromagnetic  materials  in  response  to  an  applied 
magnetic  field  can  be  explained  by  considering  two  related  mechanisms.  The  first  mechanism  is  that  domain  walls 
(the  transition  layers  between  highly  aligned  regions  termed  magnetic  domains)  move  under  the  influence  of  the 
magnetic  field  in  such  a  way  that  favorably  oriented  domains  grow  at  the  expense  of  unfavorably  oriented  domains. 
The  wall  thickness  is  determined  by  a  fine  balance  between  the  anisotropy,  which  tends  to  make  the  walls  thinner, 
and  the  Weiss- type  interaction  coupling  between  atomic  magnetic  moments,  which  tends  to  make  the  walls  thicker. 
The  second  mechanism  involves  the  rotation  of  magnetic  moments  within  domains  towards  the  field  direction. 

As  detailed  in  [7],  the  minimization  of  an  appropriate  thermodynamic  potential  A  with  respect  to  the  magnetization 
yields  the  expression 


He  = 


1  dA 

yol)M 


=  H  +  aM  + 


2/^c 


d\  da 

adM+xm 


(i) 


for  the  effective  magnetic  field,  where  fio  is  the  permeability  of  free  space.  Note  that  He  is  composed  of:  (i)  the 
applied  magnetic  field  H ,  (ii)  an  interaction  field  aM  responsible  for  the  alignment  of  neighboring  magnetic  moments, 
and  (iii)  a  component  due  to  magnetoelastic  interactions,  Ha  =  +  Aj^]. 

Maxwell-Boltzmann  statistics  can  be  used  to  quantify  the  anhysteretic  magnetization  through  the  Langevin  equa¬ 
tion 


Man  =  M, 


r  {  hb\ 

a 

coth  — 

[  V  «  J 

~  ~He_ 

(2) 


Here,  Ms  represents  the  saturation  magnetization  while  the  constant  a,  representing  the  effective  domain  density,  is 
treated  as  a  parameter  to  be  identified.  Although  the  Langevin  equation  works  well  for  paramagnetic  materials,  it  does 
not  provide  a  good  description  of  the  magnetization  phenomenon  in  ferromagnetic  materials,  in  which  the  motion  of 
domain  walls  is  impeded  by  imperfections.  These  imperfections  or  pinning  sites,  usually  attributed  to  crystal  defects 
and  the  presence  of  inclusions  such  as  dysprosium  in  Terfenol-D,  form  energy  wells  which  are  energetically  favorable 
for  domain  wall  attachment. 

The  effects  of  pinning  sites  on  domain  wall  motion  under  constant  stress  conditions  have  been  assessed  in  [8,  9,  10], 
through  consideration  of  reversible  Mrev  and  irreversible  M,rr  components  of  the  magnetization.  For  low  magnetic 
field  intensities  about  an  equilibrium  level,  the  domain  walls  bend  reversibly  while  remaining  attached  to  pinning 
sites.  As  the  applied  magnetic  field  is  increased,  the  domain  walls  achieve  sufficient  energy  to  break  free  from 
pinning  sites  while  moving  up  the  energy  well  where  they  were  originally  located,  and  attach  to  remote  sites  where 
the  energy  configuration  is  favorable.  The  effect  of  the  energy  losses  caused  by  pinning  of  domain  walls  manifests 
itself  as  hysteresis  in  the  magnetization  process. 
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Energy  balancing  yields  an  expression  for  the  rate  of  change  of  irreversible  magnetization  with  applied  field  [13], 


OMirr  Man  Mirr  ( 0He  \ 
dH  =  l ' ~dH  ) 

l—c  x  y 


(3) 


where  the  constant,  k  quantifies  the  energy  needed  to  break  pinning  sites,  c  is  a  reversibility  coefficient  which  quantifies 
the  flexibility  of  the  domain  walls,  and  S  has  the  value  +1  when  ^  >  0  and  -1  when  <  0.  The  parameter  S 
ensures  that  the  energy  lost  to  pinning  always  opposes  changes  in  the  magnetization.  Differentiation  of  (1)  with 
respect  to  H ,  along  with  application  of  the  chain  rule  yields 


dHe 

dH 
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To  facilitate  model  implementation,  the  magnetic  coupling  and  magnetoelastic  interactions  are  combined  in  a 
single  unitless  term  which  needs  to  be  identified  upon  characterization  of  A  —  M . 


3  d  (  d\  \  3  dcr  d\  3  d  (  dcr  \ 

«(M,  <r)  =  a  +  -a-  j  +  ~m~dM  +  )  • 


(5) 


The  substitution  of  (5)  into  (4), 
dM  ^  1,  yields  the  relation 


and  (4)  into  (3),  application  of  the  chain  rule,  and  recognition  from  [13]  that 


dM 


dMj  rr  Man  Mi  rr 

dH  ~  ^  -  5(M,  a) (Man  -  Mirr)  ’ 


(6) 


which  characterizes  the  irreversible  component  of  the  field-dependent  magnetization. 

Following  ideas  presented  in  [10],  the  reversible  component  of  the  field-dependent  magnetization  is  taken  to  be 
proportional  to  the  displacement  from  the  anhysteretic  magnetization,  which  may  be  expressed  in  differential  form 
as  follows 

d M r(r  f  Man  d.Mirr  A  . 

~dH~  =  c\'W~~air)  ,7) 


where  the  coefficient  c  can  be  estimated  from  the  ratio  of  the  initial  and  anhysteretic  susceptibilities  [14]  or  through 
fit  to  experimental  data. 

Finally,  the  reversible  and  irreversible  terms,  (7)  and  (6),  are  added  to  give 

d  .17  .  .  4/ a ,,  Mjrr  Man  . 

HH~[  ~C)'  ^-c  -  a (M,  c)  (Man  -  Mirr)  +  °  '  ~dlf:  1  ’ 

which  characterizes  the  total  magnetization  arising  from  the  application  of  a  magnetic  field. 


2.2  Magnetostriction  Model 

Since  the  magnetostriction  is  structure  sensitive,  assessing  a  general  A  —  M  relationship  is  not  practical  in  the  context 
of  the  present  transducer  model,  despite  the  fact  that  the  saturation  magnetostriction  As  can  be  explicitly  identified 
from  Am  and  Aioo  (Lee  [15]  provides  an  extensive  review  of  the  magnetostriction  of  anisotropic  materials.)  However, 
there  is  one  simple  case  in  which  the  relationship  A  —  M  can  be  identified,  namely  when  the  magnetic  moments  are 
forced  to  align  perpendicular  to  the  direction  of  applied  field.  This  is  accomplished  in  the  case  of  polycrystalline 
materials  such  as  Terfenol-D  by  proper  crystal  growth  during  the  manufacturing  process  and  by  applying  sufficiently 
large  compressive  preloading  during  operation.  In  this  case  the  dominant  magnetization  process  is  domain  rotation, 
and  energy  minimization  yields  the  single- valued  quadratic  law,  A (M)  =  | -^M2 . 

An  alternative  but  not  totally  unrelated  approach  consists  of  formulating  an  empirical  model  based  on  a  series 
expansion  symmetric  about  M  =  0, 


A (M)  =  7 iM2i 

i-0 


(9) 
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where  the  coefficients  7;  are  determined  from  experimental  data  [5] .  The  series  is  often  truncated  after  /  -  -  I .  yielding 
a  model  equivalent  to  the  quadratic  law  when  71  =  |-A%3  or  after  i  =  2,  yielding  a  quartic  model.  The  parameters 
71  and  72  of  the  quartic  model  are  identified  by  solving  simultaneously  the  constraints  at  magnetic  saturation  Ms 
and  at  the  inflection  point  M0, 


=  r-W;  i  T-i-w;'  =  as, 

=  271  +  1272M0  =  0. 

A  more  sophisticated  model  capable  of  including  hysteresis  effects  is  that  formulated  by  Sablik  and  Jiles  [16]. 
However,  for  the  remainder  of  this  discussion  the  quartic  magnetostriction  model  is  employed.  The  constant  term 
(i  =  0)  will  be  ignored  since  it  corresponds  to  the  elastic  strains  considered  in  Section  2.3. 


(a)  A  (Ms ) 

(b)  £Mm0) 


2.3  Transducer  Structural  Dynamics 

The  magnetostriction  A  given  by  (9)  quantifies  the  reorientation  of  magnetic  moments  towards  the  direction  of 
the  applied  field.  As  detailed  in  [7],  this  magnetostriction  is  analogous  to  the  term  d^H  in  linear  models  but  is 
inherently  nonlinear  and  hvsteretic  through  both  the  magnetization  M  and  A  —  M  itself.  It  ignores,  however,  the 
elastic  properties  of  the  magnetostrictive  material  as  it  vibrates,  represented  in  the  linear  models  by  sH a.  In  this 
section,  a  PDE  system  is  formulated  which  models  the  elastic  response  of  the  magnetostrictive  material  and  other 
transducer  components.  The  input  to  this  PDE  is  formulated  through  the  magnetostriction  A,  which  constitutes 
an  ‘internal  force’  commanding  the  vibrations  of  the  transducer.  The  solution  to  the  PDE  is  the  longitudinal 
displacements  u(t7x)  relative  to  the  prestressed  position.  Additional  details  regarding  this  PDE  formulation  are 
provided  in  [7]. 

The  structural  dynamics  are  modeled  through  consideration  of  the  magnetostrictive  rod  (5),  prestress  bolt  (6), 
spring  (7),  and  mass  load  (8)  for  the  transducer  depicted  in  Figure  1.  The  prestress  bolt  provides  a  stress  cr0  <  0  by 
compressing  the  magnetostrictive  core  against  the  spring,  modeled  by  a  linear  spring  ki  and  dashpot  c 1.  The  rod 
is  assumed  to  have  length  L,  cross  sectional  area  A,  and  longitudinal  coordinate  x.  The  material  density  is  p,  the 
elastic  modulus  is  E ,  and  the  internal  damping  is  cjj.1  Finally,  the  external  load  is  modeled  by  a  point  mass  in /, . 

The  total  stress  in  the  rod  is  given  by  the  expression 

du  32  u 

a(t,  x)  =  E— (t,  x)  +  x)  -  EX(t,  x)  +  a0  (10) 


where  the  terms  on  the  right  hand  side  represent  respectively  the  linear  elasticity  at  small  displacements,  Kelvin- 
Voigt  damping,  magnetostriction-derived  stress,  and  prestress.  Force  balancing  then  yields  the  dynamic  model  for 
the  longitudinal  displacements  and  strains. 

For  implementation  purposes,  the  model  is  formulated  in  weak  or  variational  form  by  multiplying  the  strong 
form  by  test  functions  <f>  followed  by  integration  throughout  the  length  of  the  rod.  This  reduces  the  smoothness 
requirements  on  the  finite  element  basis  since  displacements  and  test  functions  need  be  differentiated  only  once 
compared  to  the  second  derivatives  present  in  the  strong  form.  The  space  of  test  functions  is  V  =  77^(0,  L)  =  {<£>  6 
H1( 0,L)  |<p(0)  =  0},  so  that  for  all  4>{x)  6  V, 


f 


V  U 

pA-^-(t,  x)<fr(x)dx 


0u  u  u 

EA  —  (t,  x)  +  cDA-^-^{t,  x)  -  EAX(t,  ®) 


dii.  d2u 

khuyt ,  L)  +  cL  —  (t,  L)  +  L) 


m- 


(n) 


The  solution  u(t,x )  to  (11)  defines  the  longitudinal  displacements  about  the  prestressed  position.  Once  the 
displacements  are  computed,  the  strains  are  evaluated  by  taking  derivatives  with  respect  to  position,  e(t,  x)  =  |^-(t,  x), 
and  the  material  stresses  cr(t,x)  are  calculated  directly  from  (10).  Note  that  the  stress  at  the  rod  end  a(t,L)  may 

be  equivalently  calculated  from  the  boundary  condition  cr{t,  L)  =  -j  I"—  kiu{t,  L)  —  cj,  §j(0  L)  —  mL  ■§/t(^  £)j . 


1  The  parameter  E  lies  between  the  elastic  modulus  at  constant  H ,  ,  and  at  constant  B ,  EB .  Since  EB  and  EB  depend  upon  the 

field  intensity  [17],  so  does  E.  However,  for  simplicity  it  is  assumed  that  E  represents  a  nominal  or  operational  material  stiffness. 
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3  Magnetomechanical  Effect 

In  this  section,  we  consider  the  stress-dependent  magnetization  through  its  corresponding  derivative  A  uni¬ 
fying  description  of  the  changes  in  magnetization  due  to  the  action  of  stress  has  been  developed  recently  in  [4,  5,  18]. 
In  the  model  theory  developed  in  [5]  and  implemented  here,  the  main  mechanism  responsible  for  the  magnetome¬ 
chanical  effect  is  the  unpinning  of  domain  walls  that  is  caused  by  the  application  of  stress.  This  unpinning  allows 
the  walls  to  move  creating  a  corresponding  magnetization  change.  This  magnetization  change  is  such  that  the  total 
magnetization  approaches  the  anhysteretic  condition,  in  which  the  domain  configuration  has  achieved  a  state  of 
energy  minimization.  This  effect  is  depicted  in  the  schematic  M  —  H  loop  of  Figure  3(a),  where  the  magnetization 
in  the  ferromagnetic  sample  is  measured  for  varying  compressive  stresses  (while  leaving  the  magnetic  field  constant) 
starting  at  locations  above  and  below  the  anhysteretic,  points  A  and  B  respectively. 

A  second,  related  effect  is  the  change  in  the  anhysteretic  curve  itself  under  the  action  of  stress.  Application  of 
tensile  stress  produces  an  increase  in  both  the  slope  of  the  M  —  H  loop  and  the  remanent  magnetization  value, 
while  compressive  stress  produces  a  shearing  of  the  M  —  H  loop.  This  effect  is  illustrated  in  Figure  4,  where  the 
influence  of  stress  on  both  the  anhysteretic  and  total  magnetization  is  considered  for  four  stress  levels:  -10,  0,  +10, 
and  +20  MPa. 

The  net  result  of  both  effects  considered  in  combination  is  that  a  monotonically  increasing  compressive  stress  causes 
the  value  of  X  to  be  approached  (Figure  3(a)),  while  the  value  of  X  itself  decreases  (Figure  4).  The  magnetomechanical 
effect  model  presented  in  this  section  is  then  formulated  in  the  context  of:  (i)  the  effect  of  stress  on  the  anhysteretic 
magnetization  and  (ii)  a  law  of  approach  to  the  anhysteretic  magnetization  upon  application  of  stress. 


H  [A/m]  xio4 

(a)  (b) 


Figure  3:  Schematic  representation  of  the  approach  to  the  anhysteretic  magnetization  under  stress  (after  Pitman  [4]).  (a) 
Characteristic  M  —  H  loop  showing  the  approach  from  positions  above  and  below  the  anhysteretic  (A  and  B,  respectively) 
for  fixed  H.  In  either  case,  M  moves  towards  point  X  on  the  anhysteretic  curve,  (b)  Detail  of  the  variation  in  magnetization 
AM/MS  under  application  of  stress  starting  at  point  A  above  the  anhysteretic  magnetization  (A  — »  X),  and  at  point  B  below 
the  anhysteretic  magnetization  (B  — >  X). 

3.1  Anhysteretic  Magnetization 

The  anhysteretic  magnetization  calculated  in  (2)  should  be  interpreted  as  a  function  of  both  H  and  the  particular 
value  of  the  prevailing  magnetization  M ,  with  both  operating  through  the  effective  field  He  given  by  (1).  For 
instance,  at  point  A  in  Figure  5,  the  material  is  in  a  state  of  equilibrium  in  the  presence  of  zero  applied  field,  due  to 
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H  [A/m]  xio4  H  [A/m] 


Figure  4:  Model  simulations  representing  the  effect  of  stress  on  the  total  magnetization  ( - ),  and  on  the  global  anhysteretic 

magnetization  (—•—),  for  +10,  0,  -10,  and  -20  MPa. 


the  coupling  between  magnetic  moments  aM  and  the  magnetoelastic  interactions  Ha  in  the  material.  Application  of 
a  small  magnetic  field  perturbs  this  equilibrium  slightly,  giving  a  new  value  of  Man .  Applying  small  field  increments 
successively,  the  two-valued  local  anhysteretic  curve  shown  in  Figure  5  is  obtained.  Note,  however,  that  global 
equilibrium  is  never  achieved  in  this  case,  because  coupling  and  the  presence  of  pinning  sites  do  not  allow  sufficient 
relaxation. 

In  contrast,  the  ‘global’  anhysteretic  magnetization  represents  the  locus  of  points  in  which,  for  a  given  field  level 
H  and  elastic  state  a,  the  material  reaches  an  equilibrium  associated  with  the  minimization  of  the  total  energy.  The 
global  anhysteretic,  M an ,  is  a  single-valued  function  of  H  as  illustrated  in  Figure  5,  and  is  the  foundation  of  the 
magnetomechanical  effect  model  implemented  here.  It  is  noted  that  in  the  absence  of  magnetic  coupling  (cv  =  0), 
the  local  and  global  anhysteretic  curves  coincide. 

The  global  anhysteretic  depends  significantly  on  the  stress  state  in  the  material.  This  dependency  is  posed  in  terms 
of  the  effective  field  He,  by  recognizing  from  [5]  that  the  global  anhysteretic  magnetization  under  a  field  H  and  stress 
a  is  equivalent  to  the  global  anhysteretic  magnetization  under  field  He  and  zero  stress,  Man(H,a)  =  Man(He,  0). 
Assuming  the  quartic  magnetostriction  model  and  substituting  into  (1),  yields  the  effective  field 

He  =  H  +  aM  + 

2no 

Upon  definition  of  a  new  effective  coupling  term 

(27i  +  Mj-,M-)(r  +  (71  +  72.\/2).\7 

the  effective  field  may  be  written  in  the  form 


a(M,  a)  =  a  +  - — 
2po 


(27  .1/  +  472M3)<t  +  (7lM2  +  72M4) 


da 


He  =  H  +  a{M ,  a)M. 

The  global  anhysteretic  is  then  the  locus  of  points  which  satisfy  identically  the  modified  Langevin  equation  (2), 

' H  +  a(M an,  a)M a 


coth 


H  +  a(M an,  a)M a 
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H  [A/m]  xio4 

Figure  5:  Model  simulations  showing  the  relative  magnetization  M/Ms,  relative  local  anhysteretic  Man(H,  M,  a)/Ms  = 
coth[(i/  a(M,  a)M)/a]  —  [a/(H  -)-  a(M,  cr)M)],  and  relative  global  anhysteretic  Man{H,a)/Ms  =  coth[(if 

o(Mo»,  v)Man)/a\  —  [a/(H  +  ct(Man,cr)Man)],  as  a  function  of  H. 

3.2  Stress-Dependent  Magnetization 

The  behavior  of  the  magnetostrictive  core  depends  on  more  than  just  the  externally  applied  stress  a  and  magnetic 
field  H .  In  fact,  the  magnetic  history  of  the  specimen,  or  equivalently  the  displacement  from  the  anhysteretic 
condition  in  the  case  of  major  M  —  H  loops,  is  what  determines  the  direction  of  changes  in  magnetization.  According 
to  the  experimental  observations  documented  in  [3,  8,  19],  the  direction  in  which  the  magnetization  changes  with 
applied  stress  is  independent  of  the  sign  of  the  stress,  for  small  stresses  and  when  the  magnetization  is  sufficiently 
distant  from  the  anhysteretic.  It  is  then  inferred  that  the  direction  of  change  is  dependent  not  on  the  stress  itself, 
but  on  a  quantity  which  is  independent  of  the  sign  of  the  stress.  In  this  context,  Jiles  hypothesized  in  [5]  that  this 
quantity  is  the  elastic  energy  per  unit  volume,  W  =  a2/{2E),  which  is  clearly  independent  of  the  sign  of  a  and  which 
may  be  responsible  for  some  unpinning  of  domain  walls.  The  ‘law  of  approach’  to  the  anhysteretic  condition  is  then 
formulated  as  follows:  the  rate  of  change  of  magnetization  with  elastic  energy  is  proportional  to  the  displacement  of 
the  prevailing  magnetization  from  the  anhysteretic  magnetization.  As  before  (Section  2.1),  this  may  be  posed  via 
irreversible  and  reversible  components.  The  law  of  approach  as  applied  to  the  irreversible  component  then  gives 

=  \(M„,  -  */.„),  (12) 

where  £  is  a  coefficient  to  be  identified  with  dimensions  of  energy  per  unit  volume.  Making  the  substitution  4^-  =  ^ 
in  (12),  along  with  application  of  the  chain  rule,  permits  writing  the  rate  of  change  of  Mirr  with  cr  as  follows, 

dMirr  <7  , - 

Qa  =  —  Mirr).  (13) 

A  procedure  similar  to  that  employed  in  the  field-dependent  case  yields  the  rate  of  change  of  reversible  magneti¬ 
zation  with  stress, 


BMrev  f  9M an  0  Mjrr  \ 

da  \  da  da  )  ’ 


(14) 


in  which  c  is  the  coefficient  defined  in  (7). 


The  reversible  and  irreversible  terms,  (14)  and  (13),  are  added,  to  give 


dM 

da 


(1  -c)~(Man-MirT)  +  c^^. 


which  characterizes  the  magnetization  arising  from  the  application  of  stress. 


(15) 


4  Coupled  Model 

In  the  presence  of  both  an  applied  magnetic  held  H  and  stress  a:  the  total  magnetization  is  dictated  by  the 
superposition  of  the  held-  and  stress-dependent  components  given  in  (8)  and  (15).  For  implementation  purposes, 
however,  we  consider  instead  the  governing  differential  equation 

(t Mjrr  /  0  M ;r  r  \  d  1 /  /  dMirr  \  da 

dt  \  dH  J  dt  \  da  )  dt 

for  the  total  irreversible  magnetization  Mjrr  which  arises  from  the  superposition  of  held  and  stress. 

It  is  noted  that  the  relationship  between  applied  current  I(t )  and  magnetic  held  H{t,  x )  generated  in  a  solenoid  is 
often  characterized  by  H  (t)  =  nl(t),  where  n  =  Nl°e^g[^'s  ■  Experimental  evidence  on  research  transducers  indicates 
that  H  =  nl  does  not  provide  a  sufficiently  accurate  characterization  of  the  magnetic  held,  so  an  empirical  position- 
dependent.  filter  <3?(a,’)  is  considered,  which  accounts  for  end  effects,  demagnetizing  factors,  and  overall  losses  in  the 
magnetic  circuit.  Thus,  the  time  and  spatial  dependencies  of  the  held  are  formulated  via  H(t,x)  =  /(!)$(*). 2 
Applying  the  chain  rule  and  substituting  (6)  and  (12)  into  (16)  yields 


dMirr 

dt 


{t,x) 


Man  (t,  x)  Mirr{t:X^j 


-  a(M,  a)(Man(t ,  x)  -  Mirr(t,  x)) 


dl 

dt 


W  •  *(*) 


a  (t,  x)  ,  ,  . ,  ,  ..da , 

'  ( M an(t ,  x)  —  Mirr{t,  a.’))  •  — (E  *), 


(17) 


which  after  integration  gives  the  irreversible  component  of  the  magnetization,  M,rr(f,  x). 

The  reversible  magnetization  is  calculated  directly  upon  integration  and  subsequent  superposition  of  (7)  and  (13), 

Mrev  (t r  a)  —  C  [Afan(f,  x)  T  M  an(t ,  a)J  cMjrr  (f ,  x) . 


The  total  magnetization  arising  from  application  of  both  a  held  and  stress  is  then 

M(t ,  x)  =  c  [Man  (t,  x)  +  Man(t ,  x)]  +  (1  -  c)Mirr{t ,  x), 

which  incorporates  the  reversible  and  irreversible  components  of  magnetization.  It  should  be  noted  that  in  the 
case  of  constant  stress  (^j  =  0)  or  constant  held  (-|y  =  0),  the  expression  reduces  to  the  individual  components 
characterized  in  (8)  and  (15). 


5  Model  Validation 


5.1  Experimental  Setup 

In  this  section,  the  model  presented  in  Section  4  is  applied  to  a  transducer  developed  at  Iowa  State  University 
(the  configuration  of  which  follows  the  general  guidelines  of  Figure  1),  and  used  to  characterize  the  time-dependent 
magnetization  and  displacement  produced  by  the  transducer  in  response  to  an  applied  current  I(t).  These  calculations 
are  compared  to  experimental  measurements  collected  in  this  transducer.  The  driver  was  a  4.53  in  (115  mm)  long, 
0.5  in  (12.7  mm)  diameter,  monolithic  Tb0 .3DV0 ^Fei  9  rod  manufactured  using  the  Bridgman  process.  Surrounding 

2Note:  (i)  't ( .?;)  is  equal  to  n  in  the  idealized  situation  of  an  infinitely  long,  lossless  solenoid,  and  (ii)  an  alternative  approach  for 
identification  of  H  —  I  consists  of  solving  the  Ampere’s  law  or  the  Biot  Savart  law.  Several  software  packages  utilize  finite  element  analysis 
to  accomplish  this. 
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the  rod  were  an  innermost  single-layer,  150-turn  pickup  coil,  used  for  magnetic  induction  measurements,  and  a  six- 
layer,  900-t.urn  drive  coil  connected  to  a  Techron  7780  current  control  amplifier.  Both  solenoids  were  wound  in-house 
with  26  AWG  magnet  wire. 

Steel  end  caps  and  a  cylindrical  Alnico  V  permanent  magnet  provided  flux  closure.  For  the  experiments  reported 
here,  the  permanent  magnet  was  demagnetized  to  obtain  unbiased  operation.  The  displacement  plunger  was  guided 
by  a  nylon  bushing  mounted  within  an  aluminum  head,  which  was  held  in  place  by  noil-magnetic  steel  bolts  attached 
to  a  heavy  stainless  steel  base.  The  0.7  Hz  sinusoidal  reference  signal  input  to  the  amplifier  was  provided  by  a 
Tektronix  2642A  Fourier  Personal  Analyzer,  which  was  also  used  to  record  the  pick-up  coil,  drive  current,  drive 
voltage  and  displacement  signals.  The  external  load  was  a  0.5  kg  mass  threaded  onto  the  displacement  plunger. 

The  compressive  preloading  necessary  to  avoid  tensile  stresses  on  the  moderately  brittle  Terfenol-D  rod  was 
provided  by  a  steel  bolt  located  within  the  transducer  base.  The  bolt  pushed  the  rod  against  a  Belleville  washer,  whose 
force-displacement,  characteristic  curve  was  measured  in  an  MTI  Phoenix  compression  machine.  During  preloading, 
the  bolt  was  tightened  until  the  plunger’s  displacement  matched  the  displacement  in  the  washer’s  characteristic  curve 
corresponding  to  the  desired  force  level. 

The  displacement  at  the  rod  tip  was  measured  with  a  Lucas  Schaevitz  LVM-110  linear  variable  differential  trans¬ 
ducer  equipped  with  a  050-MHR  probe,  capable  of  99.9923  %  linearity  over  ±0.050  in.  The  corresponding  bulk 
strain,  used  for  comparison  to  model  simulations,  was  computed  by  dividing  this  displacement  by  L. 

For  identification  of  the  $(*)  functional,  a  series  of  flux  measurements  was  conducted  with  a  Hall  effect  probe 
(F.  W.  Bell  9500  series),  located  within  the  transducer  to  allow  characterization  of  end  effects,  demagnetizing  factors, 
and  magnetic  circuit  nonlinearities.  The  amplitude  of  the  flux  measured  adjacent  to  the  surface  of  the  rod  at  locations 
3  mm  apart  over  a  longitudinal  line,  \Bsur(t,x)\,  is  equal  to  the  magnetic  field  amplitude,  \Hsur(t,  a?) | ,  because  in 
air  B  =  H  in  the  CGS  units  system.  The  filter  was  identified  by  taking  the  ratio  at  several  current  levels, 

and  a  look-up  table  was  used  at  the  desired  current  level. 

The  magnetic  induction  was  calculated  by  integration  of  the  pickup  coil  signal  Vpu,  following  the  Faraday-Lenz 
law  of  magnetic  induction,  B(t)  =  —  N  lA  f*  VPu(r)dT.  The  parameter  Apu  is  the  mean  cross  sectional  area  of  the 
pickup  coil,  and  Npu  was  equal  to  150  turns  as  indicated  before.  The  magnetization  used  for  comparison  against 
model  simulations  was  then  computed  from  the  magnetic  constitutive  equation  M  =  j^-B  —  H . 

5.2  Parameter  Estimation 

A  total  of  fourteen  parameters  completely  describe  the  coupled  magnetic-structural  model,  of  which  six  are  magne¬ 
tization  parameters  (a:k,c,a,  Ms  and  £;),  two  are  magnetostriction  parameters  (M0  and  As,  or  71  and  72),  three  are 
material  properties  of  the  transducer  driver  ( E ,  p  and  cp),  and  three  describe  the  load  characteristics  (&l,C£  and 
rriL).  The  parameters  were  quantified  in  three  manners:  i)  by  direct  measurement  (&l,  niL  and  p),  ii)  from  published 
specifications  for  Terfenol-D  ( Ms ,  E  and  As),  and  iii)  by  optimization  through  a  least  square  fit  to  experimental  data 
(a,  k,  c,  a,  Mo,  cd  ,  cl  and  !;).  The  optimization  was  performed  in  two  steps.  In  the  first,  the  five  magnetization  pa¬ 
rameters  and  damping  components  were  estimated  from  magnetization  data.  The  second  step  consists  of  identifying 
the  magnetostriction  coefficients  from  strain  data. 

The  magnetization  parameters  have  been  previously  identified  for  ferromagnetic  materials  in  general  in  [14],  and 
for  Terfenol-D  in  particular  in  the  uncoupled  case,  in  [6,  20].  The  parameter  optimization  performed  for  the  coupled 
case  reveals  some  differences  in  these  four  parameters  compared  to  the  previous  work,  which  is  most  likely  attributable 
to  the  effects  of  coupling. 

For  Terfenol-D,  the  nominal  saturation  magnetostriction  is  As  =  1000  x  10“ b.  The  saturation  magnetization  has 
been  found  to  deviate  little  from  the  nominal  Ms  =  0.765  x  10®  A/m.  The  optimized  model  parameters  employed  in 
the  simulations  are  summarized  in  Table  1. 

5.3  Magnetization  and  Strain  Results 

The  mean-field  model  discussed  in  Section  2.1  provides  a  characterization  of  the  magnetization  M  in  the  transducer 
core  in  response  to  an  applied  field  H .  The  performance  of  the  model  is  demonstrated  in  Figure  6(a).  It  is  noted 
that  the  model  accurately  characterizes  the  fundamental  aspects  of  the  magnetization  characteristic  curve.  However, 
the  model  does  not  quantify  the  constricted  behavior  at  low  field  levels,  which  is  perhaps  due  to  domain  rotation 
effects.  It  is  observed  that  the  magnetization  model  does  not  include  the  moderate  anisotropy  typical  of  Terfenol-D. 

The  strain  simulations  are  gauged  against  experimental  data  in  Figure  6(b).  The  model  accurately  captures  the 
nonlinearity,  hysteresis,  and  saturation  features  present  in  the  data.  We  note  that  the  errors  observed  in  the  M  —  H 


10 


model  propagate  to  the  strain  measurements  through  the  magnetostriction  law,  a  fact  which  explains  some  of  the 
discrepancies  present  in  the  performance  of  the  complete  e  —  H  model.  Moreover,  the  approximate  nature  of  the 
quartic  law  itself  is  thought  to  constitute  an  additional  source  for  errors. 


(a)  (b) 

Figure  6:  Comparison  of  model  simulations  ( - )  with  experimental  data  (—•—),  under  a  prestress  of  1  ksi  (6.9  MPa),  (a) 

Magnetization  of  the  magnetostrictive  core,  (b)  total  strain  output  by  transducer. 


Magnetization 

Magnetostriction 

Structural 

Load 

Geometry 

a  =  6500  A/m 
k  =  6000  A/m 
c  =  0.008 
a  =  0.03 

Ms  =  0.765  x  106  A/m 
£  =  24.5  x  103  Pa 

M0  =  0.8  x  10l:i  A/m 

As  =  1000  x  10-6 

71  =  2.9  x  10“ 15  nr/A2 

72  =  -7.5  x  10-28  m4/A4 

E  =  30  x  10y  Pa 
p  =  9250  kg  /  m3 
cj  =1  x  106  Ns/rn 

kL  =  2.5  x  10“  N/m 

niL  =  0.5  kg 

cl  =  lx  103  Ns/rn 

L  =  0.115  m 

D  =  0.0127  m 

Table  1:  Physical  parameters  and  dimensions  used  in  the  coupled  elastic-magnetic  model. 


6  Concluding  Remarks 

A  coupled  elastic-magnetic  model  for  the  magnetization  and  strain  response  of  magnetostrictive  transducers  in 
response  to  applied  currents  has  been  presented.  The  model  simulations  accurately  predict  major  loop  behavior  under 
mechanical  preloads  and  quasi-static  frequencies  of  operation.  As  indicated  before  [6],  the  model  also  accurately 
models  the  transducer’s  fundamental  frequency  within  5%  of  its  measured  value.  Nevertheless,  the  model  in  its 
present  format  should  be  restricted  to  near-DC  frequencies.  Work  is  in  progress  to  include  AC  losses  and  to  extend 
the  model  to  the  frequency  domain. 

It  is  desirable  to  take  advantage  of  the  force  simulations  that  the  coupled  model  is  able  to  generate.  Work  is 
in  progress  to  gather  experimental  information  for  validation  of  the  forces  predicted  by  the  model.  Preliminary 
indications  suggest  that  the  involved  stress  amplitudes  at  input  levels  similar  to  those  utilized  in  Section  5  are  on  the 
order  of  3.2  MPa,  or  two  orders  of  magnitude  below  the  stress  levels  used  in  the  experimental  works  of  Pitman  [4]  and 
the  theoretical  developments  of  Jiles  [5].  This  fact  explains  why,  for  the  operating  regimes  typical  of  the  particular 
transducer  unit  tested,  the  coupled  model  did  not  provide  a  significant  improvement  over  the  uncoupled  model 
presented  previously  [6]. 

We  believe  that  the  added  accuracy  and  flexibility  of  the  elastic-magnetic  coupled  model  will  enhance  the  control 
and  performance  capabilities  of  smart  structure  systems  involving  magnetostrictive  materials,  in  particular  those 
utilizing  the  full  compression  range  of  Terfenol-D,  of  about  700  MPa.  The  uncoupled  formulation  will  still  prove 
highly  useful  in  applications  where  stress  variations  are  less  significant. 
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